A&A manuscript no. 

(will be inserted by hand later ) 



Your thesaurus codes are: 

06(02.01.1; 02.18.5; 02.19.2; 08.14.1; 08.16.6) 



ASTRONOMY 

AND 
ASTROPHYSICS 



Electron energy losses near pulsar polar caps: 
a Monte Carlo approach 



o 
o 
o 

(N 

c 

o\ 

(N 
> 

m 

On 
O 
CO 

o 

O 
O 

33 

Oh' 

6 

CO 



X 



J. Dyks and B. Rudak 

Nicolaus Copernicus Astronomical Center, Rabiariska 8, 87-100 Toruri, Poland 
E-mail: jinx@ncac.torun.pl, bronek@camk.edu.pl 

Received / accepted 



Abstract. We use Monte Carlo approach to study the en- 
ergetics of electrons accelerated in a pulsar polar gap. As 
energy-loss mechanisms we consider magnetic Compton 
scattering of thermal X-ray photons and curvature radia- 
tion. The results are compared with previous calculations 
which assumed that changes of electron energy occurred 
smoothly according to approximations for the average en- 
ergy loss rate due to the Compton scattering. 

We confirm a general dependence of efficiency of elec- 
tron energy losses due to inverse Compton mechanism on 
the temperature and size of a thermal polar cap and on the 
pulsar magnetic field. However, we show that trajectories 
of electrons in energy-altitude space as calculated in the 
smooth way do not always coincide with averaged Monte 
Carlo behaviour. In particular, for pulsars with high mag- 
netic field strength (B pc ;> 3 x 10 12 G) and low thermal 
polar cap temperatures (T <J 5 x 10 6 K) final electron 
Lorentz factors computed with the two methods may dif- 
fer by a few orders of magnitude. We discuss consequences 
for particular objects with identified thermal X-ray spec- 
tral components like Geminga, Vela, and PSR B1055— 52. 

Key words: acceleration of particles - radiation mecha- 
nisms: non-thermal - scattering - stars: neutron - pulsars: 
general 



1. Introduction 

Curvature radiation (CR) and magnetic inverse Compton 
scattering (ICS) are usually considered to be the most 
natural ways of hard gamma-rays production operating 
at the expense of pulsar rotational energy (eg. Zhang & 
Harding 2000, and references therein). These two radia- 
tion mechanisms dominate within two different ranges of 
Lorentz factors 7 of beam particles (ie. those leaving the 
polar cap). When 7 < 10 6 , magnetic inverse Compton 
scattering plays a dominant role in braking beam parti- 



cles (Xia et al. 1985; Chang 1995; Sturner 1995) and is 



When 7 ;> 10 6 , the curvature radiation becomes re- 
sponsible for cooling beam particles (eg. Daughcrty & 
Harding 1982), yet, the inverse Compton scattering may 
still be important for secondary e ± -pair plasma. Zhang & 
Harding ( |2000| ) show that different dependencies of X-ray 
and gamma-ray luminosities on the pulsar spin-down lu- 
minosity as inferred from observations by CGRO, ROSAT 
an d ASC A (Arons |1996t Becker & Triimper [I997t Saito et 



al. 1998, respectively) can be well reproduced when this 



effect is included. 

So far, the influence of ICS on electron energ y hav e 
been inv estiga ted most thoroughly b y Cha ng ( 1995 ), 
Sturner fll995| ), Hardi ng & Muslimov ( |l998|) , and Sup- 
per & Triimper (200C). One of several interesting results 
they found was that energy losses due to resonant ICS 
can limit the Lorentz factors 7 of electrons to a value 
which depends on both electric field strength E\\, tem- 
perature T and radius i?th of thermal polar cap, and on 
the magnetic field strength B pc at a polar cap. For exam- 
ple Sturner (1995|) shows that assuming the acceleration 



model of Michel ( |1974 ) and a thermal polar cap size com- 
parable to that determined by the open field lines, the 
Lorentz factors 7 are limited to ~ 10 3 if B pc > 10 13 G, 
and T > 3 x 10 6 K. This acceleration stopping effect be- 
comes more efficient for stronger magnetic fields, and it 



led Sturner (1995) to propose it as a possible explanation 



for an apparent cutoff around 10 MeV in gamma radiation 
of B1509-58 (Kuiper et al. 



1999) 



the mai n sou rce of hard ga mma- ray photons (Sturner & 
Dermer |l994| ; Sturner et al. |l995[ ) . 



Those results were obtained with a numerical method 
which assumes smooth changes in electron energy and 
makes use of expressions fo r ave raged energy losses due to 
the magnetic ICS (Dermer 199C). The continuous changes 
of the Lorentz factor with height as determined with this 
fully deterministic treatment are considered as representa- 
tive for a behaviour of a large number of electrons. Actu- 
ally however, electrons lose their energy in discontinuous 
scattering events. Between some of these events (which oc- 
casionally can be very distant even for a short mean free 
path) the energy of an accelerated electron may increase 
considerably (even a few times) and therefore, the elec- 
tron can find itself in quite different conditions than if it 
were losing its energy continuously (the mean free path 
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depends very strongly on 7). Since the energy loss rate 
due to the resonant ICS is not a monotonic function of 
7 this effect may be essential for future behaviour of the 
particle. 

The aim of our paper is to investigate the energetics 
of electrons on individual basis with a Monte Carlo treat- 
ment. In particular, a question of the ICS-limited acceler- 
ation is adressed. In Section 2 we calculate electron energy 
losses near a neutron star using both, the Monte Carlo ap- 
proach and then the smooth method (after Dermer 1990; 
Chang 1995; Sturncr 1995) for comparison. We consider 
two distinct cases of electron acceleration: 1) in the con- 
stant electric field, proposed by Michel (1974) (hereafter 
M74); 2) in the electric field elaborated by Harding & 
Muslimov Q199SD (hereafter HM98). Section 3 presents the 
results for electron energetics, and the finding that within 
some range of pulsars parameters, the averaged Monte 
Carlo behaviour of electrons does not coincide with a so- 
lution found with the smooth method. It also contains the 
explanation of these differences as well as consequences for 
predicted efficiencies of gamma-ray emission. Conclusions 
are given in Sect. 4. 

2. The model 

We consider electrons accelerated in a longitudinal electric 
field induced rotationally within the region adjacent to the 
surface of a neutron star. The particles lose their energy 
due to scattering off soft thermal X-ray photons through 
magnetic inverse Compton mechanism and (marginally) 
due to emission of curvature radiation. The thermal pho- 
tons are assumed to originate from a flat 'thermal polar 
cap' with a temperature T ~ a few x 10 6 K and with a 
radius R ta — a few x r pc where r pc = (27r R^/cP) 1 ^ 2 is 
the standard polar cap radius as determined for an aligned 
rotator by the open lines of purely dipolar magnetic field; 
(R ns denotes the neutron star radius and P is the pulsar 
period). 

Since the ICS losses are to dominate over the curva- 
ture losses, one has to ensure that electrons do not attain 
extremely high Lorentz factors (7 < 10 6 ). Therefore, ci- 
ther the electric field should be weak or at least the size 
of the accelerating region should be small enough to pre- 
vent 7 ;> 10 6 . For this reason most papers focusing on the 
resonant ICS made use of a relatively weak electric field 
after Michel 1974| (for example: Sturner 1995 and recently 
Supper & Triimper |2000 ). Michel's model is considered as 
unrealistic since it ignores magnetic field line curvature 
(detailed treatment of th is pro blem was first presented 
by Arons & Scharlemann 1979 ) as well as inertial frame 
dragging (its importance was acknowledged and the effect 
was worked out for the first time by Muslimov & Tsygan 
1992| ). Both effects lead to significantly stronger electric 



fields operating in polar-cap accelerators. 

However, for the sake of comparison with previous pa- 
pers on the resonant ICS we first consider the weak electric 



field after Michel 1974. Accordingly, we assume that elec- 
trons are accelerated in a constant En extending between 
stellar surface (h = 0) and altitude h acc : 



(S7rmcB vc \ 
e p J 



1/2 



with h* 



3.58 ■ 10 5 (^''"Vcm" 1 



5 ( B 12 \ 1 / 2 - 

for < h < h acc ; 
for h > h &cc . 



(1) 



r pc , and where P is in seconds, B12 is B pc in 
Teragauss. Fixed location of h acc means that the model is 
not self-consistent in the sense that the accelerating field 
will not be shorted out by a pair formation front (should 
it occur at a height below h acc — r pc ). 

As the second case for E\\ we take the advanced model 
of an electric field in the form elaborated by Harding & 



Muslimov (1998). HM98 were the first to incorporate elec- 
tron energy losses due to the ICS in a strong- E\\ accelera- 
tion model. However, this was done to introduce a high al- 
titude accelerator with negligible resonant ICS losses (an 
idea of an unstable ICS-induced pair formation fronts). 
In that work Harding and Muslimov also calculated self- 
consistent values of the acceleration height h acc using an 
accelerating electric field which takes into account the up- 
per pair formation front (eqs. (18) and (23) in HM98). 
For h acc <J 0.1 r pc the field is well approximated with 



E\\(h) 



B, 



pc 



h h a 



C 1 — CGR Rns Rn 

x [k cos x + A sin x] , 



(2) 



where f2 = 2n/P, (1 — £gr) — 0.6, K ~ 0.15, x is the angle 
between the magnetic dipole and rotation axes, h is the 
altitude, and the factor Asmx is negligible for nonorthog- 
onal rotators (see HM98 for details). For /i acc approaching 
r pc the linear dependence on h &cc in Eq. (|^) disappears. We 
have found that for 0.5r pc < /i acc < 3r pc and h < r pc /3 
the electric field given by eq.(18) in HM98 is well aproxi- 
mated with a formula 



E\\(h) 



QR n 



B„ 



C (l-e G R) 1/2 RnsRn 

x [k cos x + A sin x] ■ 



1 - 



(3) 



For purposes of comparison as well as to simplify our 
calculations, we inject electrons from the center of polar 
cap and propagate them along the straight magnetic axis 
field line. The propagation proceeds in steps of a size dh, 
which is one hundred times smaller than any of character- 
istic length scales involved in the problem, like the acceler- 
ation length scale or the mean free path for the magnetic 
ISC. In each of the distance steps the energy of an electron 
is increased by a value eE^dh. For the sake of complete- 
ness, we also include the energy loss due to the curvature 
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radiation \j CI \mc 2 dh/v where j cr is the CR cooling rate where a T is the Thomson cross section, and gi and g^ are 



given by 



given by 



7c 



2e 2 /3 4 7 4 
3mcf% UTV 



(4) 9i(u) = 



and v = c(3 is the electron velocity. Nevertheless, the CR 
is a very inefficient cooling mechanism in the presence of 
acceleration models we consider below and it could be 
neglected equally well. When computing energy losses due 
to the CR we assume artificially (after Sturner 1995) that 
the magnetic axis has a fixed radius of curvature p curv = 
10 7 cm, though we keep the electrons all the time directly 
over the polar cap center. 

Magnetic inverse Compton scattering has been treated 
with a Monte Carlo simulations. The framework of our nu- 
merical code is based on an approach proposed by Daugh- 
erty & Harding (1989). We improve their method of sam- 
pling the parameters of incoming photon and account for 
the Klein-Nishina regime in an approximate way. 

At each step in the electron trajectory, the optical 
depth for magnetic ICS is calculated as dr — dhlZ/c, 
where 1Z is a scattering rate, to decide whether a scatter- 
ing event is to occur. The scattering rate 1Z in the observer 
frame OF is calculated as 



* = c /^/^(Sf0 (1 



(5) 



(eg. Ho & Epstein 1989) where = d/j,d(j> is the solid an- 
gle subtended by the source of soft photons, (fi = cos 8), 
a is a total cross section (see below), and dn p h/ds/d£l is 
the density of the soft photons per unit energy and per 
unit solid angle. The symbol £ denotes photon energy in 
dimensional energy units. Hereafter we will use its dimen- 
sionless counterpart e = s/(rnc 2 ) to denote the photon 
energy in the observer frame OF and the primed symbol 
e' = ej(l — j3fi) in the electron rest frame ERF. 

For temperatures of the thermal polar cap and Lorentz 
factors of electrons considered below, photon energies e' 
may fall well above e B = B/B c[ , a local magnetic field 
strength in units of the critical magnetic field B cr = 
m 2 c 3 (efr)~ 1 . This suggests a full form of the relativistic 
cross section for Compton scattering in stron g magnetic 
fields to be used (Daugherty & Harding 1986 ). However, 
incoming photons that propagate in the OF at an an- 
gle 9 = cos -1 fj, are strongly collimated in the ERF with 
a cosine of polar angle // = (/x — /3)/(l — j3fi) c lose to 
T. As was shown by Daugherty & Harding ( 1986| ), when 



approaches 1, resonances at higher harmonics become 
narrower and weaker, and scattering into higher Landau 
states becomes less important. In such conditions, the 
polarization-averaged relativistic magnetic cross section in 
the Thomson regime is reasonably well approximated with 
a nonrelativistic, classical limit: 



l-^ + d + p' 2 ) 



Hi 



9 2 - .91 



(G) 



(u + 1) 2 ' 



92 (u) = 



(u - l) 2 + a 2 



(7) 



with u = e'/e B , a = 2afe B /3, wher e a f i s a fine-structure 
constant (eg. Herold 1979| ; Dermer 1990| ). 

In the Klein-Nishina regime (e' ;> 1) the relativis- 
tic magnetic cross section for the ~ 1 case becomes 
better approximated with the well known Klein-Nishina 
relativistic nonmagnetic total cross section a KN (Daugh- 
erty & Harding [1986| Dermer |l990p . When e' > e B , (the 
condition fulfilled in the K-N regime since e B < 1 holds 
throughout this paper), the resonant term (52 — <?i)/2 in 
Eq. (^) becomes negligible and the nonresonant term g\ 
approaches unity which results in er « a T . Therefore, to 
approximate the Klein-Nishina decline we replace the sin- 
gle nonresonant term 171 in the square bracket of Eq. ^ 
with a KN /a T for e' > 2e B . 

To make calculations less time-consuming, we have 
used the delta-function approximation for the resonant 



part of the cross section (Dermer 1990), when calculating 
the optical depth for a scattering and the averaged energy 
loss of electron (Eq.[ll], see below). 

As the density of soft photons dnp^/ds/dft we take the 
spectral density of blackbody radiation given by 



dn 



ph\ 



dsdQ, J 



dsdVL 



e 2 dedtt 



A3 [exp(e/T) - 1] 



(8) 



and 



where T is a dimensionless temperature (T = ^2) 
A c = h/(mc) is the electron Compton wavelength. Hereby 
we neglect anisotropy of the thermal emiss ion ex pected in 
the strong magnetic field (eg Pavlov et al. 1994 ). For e <C 
e B a preferred photon propagation direction is that of the 
magnetic field because the opacity is reduced for radiation 
polarized across B. The probability of collision with an 
electron is greatly reduced for photons propagating at \x « 
1, thus, the efficiency of ICS as calculated below should 
be treated as an upper limit. 

If a scattering event occurs, we draw the energy e and 
the cosine of polar angle /x for an incoming photon from a 
distribution determined with the integrand of Eq. (^). We 
do this with the simple two-dimensional rejection method 
which is reliable though more time-consuming (cf. an ap- 
proach by Daugherty & Harding 1989). Next, we trans- 



form e and fi to the ERF values e' and //, and sample the 
cosine of polar angle of an outgoing photon /i' s using the 
differential form of the cross section (||) in the Thomson 
limit: 



da 3<7 T , , 
de'dil' = 16^ ^ s 



t ) (I-a^XI-^V 



+ \ (i + m' 2 )(i + a4 2 )G?i + s 2 ) 



(9) 
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(eg. Herold 1979), where dfl' s = d<fi' s dij,' s is an increment 
of solid angle into which outgoing photons with energy 
e' s in the ERF are directed. As in the case of total cross 
section (||), when e' > 2e B we replace the factor (g± + g 2 ) 
in ([)]) with (2<7 KN /cr T + g 2 — pi). The sampled value of 
/j,' s determines the energy £g Cat for the particular scattered 
photon with the relativistic formula: 



l + e'(l-MVs) + 



l + 2 e y s (^- / /) + e' 2 ( A1 ' s -^) : 



1/2 



(10) 



appropriate for collisions with a recoiled elect ron remain- 
ing at the ground Landau level (Herold 1979 ). Finally, a 
value of the electron's longitudinal momentum in the ERF 
is changed due to recoil from zero to (e'fi' — e' scat /j,' s )mc, 
and transformed back to the OF. 

The Monte Carlo method will be compared with a 
smoo th integr ation treatment of the magnetic ICS (Chang 
1995 ; Sturner 1995|) . In the integration method 



we assume 



the same procedure to accelerate an electron and to sub- 
tract its energy losses due to the curvature radiation as 
described above (Eq. ||) . The only difference is in account- 
ing for electron energy losses due to the magnetic inverse 
Compton scattering. These are estimated in each step (re- 
gardless the value of the optical depth for the scattering 
process) from the following formula for the mean electron 
energy loss rate: 



-Tics 



de / dfl 



de' / dtt' 



da 



de'dfl' 



(es-e) (11) 



where e s = e' s j(l + /3{ i' s ) is the scattered photon energy 



in the OF (eg. Dernier 199C). In other words, the electron 
Lorentz factor is determined as a solution of the differen- 
tial equation: 



d-y , 

-TT = V (7acc + 7ics + 7cr)- 

dh 



(12) 



At each step, as an electron moves upwards, a decrease 
in both the dipolar magnetic field strength and the solid 
angle subtended by the thermal polar cap is taken into 
account in both methods. 



3. Results 

We compare the Monte Carlo method described in Sect. 2 
with the integration method (numerical integration of 
Eq. (|l^)), for a pulsar with the polar magnetic field 
strength and rotation period as for B1509— 58 (B pc = 
15.8 x 10 12 G, P = 0.15 s). 



3.1. The case of the M74 electric field 

In this subsection the results are presented for the elec- 
tric field of Michel Ql974|) (see Eq.ll]). The assumed val- 



ues of B pc and P give Eu = 3.67 x 10 b V cm -1 and the 
standard polar cap radius r pc = 3.74 x 10 4 cm. Following 
Sturner ( |l995|) , we assumed the thermal polar cap radius 
Rth = 10 5 cm (~ 2.7 x r pc ) and h acc = r pc . Calcula- 
tions have been performed for four different temperatures 
of the thermal polar cap: Tq = 3.0, 3.5, 4.0 and 4.5, where 
Tq = T/(10 6 K). 

Each panel in Fig. [l] presents in grey the curves j(h) as 
calculated in the Monte Carlo way for one hundred indi- 
vidual electrons. The behaviour determined with Eq. ( p"2| ) 
is overplotted as the thick solid line. The inclined dot- 
ted line presents changes of electron energy if there were 
no losses. It becomes horizontal at altitude h acc = r pc = 
3.74 x 10 4 cm above which the longitudinal component of 
electric field En is assumed to be screened by a charge- 
separated plasma distribution (see Eq.[l]). This height is 
denoted by the vertical long dashed line. 

The thick solid trajectories j(h) in four panels of Fig. [l] 
represent the "averaged treatment" solutions and are in 
good agreem ent w ith the results of previous calculations 
by Sturner (|l995|) (cf his Fig. 4). Below h ~ 10 cm 
they overlap with the case with no energy losses (dot- 
ted), since the acceleration rate of electron (as given by 

7acc = ^?P C ) significantly exceeds the energy loss rate 
due to the ICS (see Fig. ||a). Energy losses due to the 
CR are negligible for any Lorentz factor accessible for an 
electron, given the assumed acceleration model of Michel 
(Eq. [|). At larger altitudes the thick trajectories illustrate 
the acceleration stopping effect noted by Chang (1995) 
and Kardashev et al. (1984): the energy that the electrons 
would have gained between altitudes h ~ 10 2 and ~ 10 4 
cm due to the electric field E\\ is transferred to the ther- 
mal photons through the resonant inverse Compton scat- 
terings. Depending on which process - either the resonant 
ICS cooling or the acceleration - ceases first, the elec- 
trons end up with energy which is either closer to eE\\ ha. cc 
(ie. the energy acquired due to the full voltage drop at no 
radiative losses, the case Tq = 3.0) or decreases towards a 
few tens xmc 2 (the cases with Tq > 3.5). 

As can be seen in Fig. [|, the Monte Carlo tracks gen- 
erally behave qualitatively in the same way: the increas- 
ing temperature of the thermal cap increases the altitude 
at which the acceleration takes over, and eventually, the 
Lorentz factors of most of electrons become limited below 
10 2 . However, there are strong differences between the av- 
erage energy of the Monte Carlo electrons and the value 
obtained with the integration method. They are especially 
pronounced for Tq = 3.5, which is the case where a max- 
imum altitude at which the acceleration can be counter- 
balanced by the electron energy losses due to the reso- 
nant ICS is equal to h acc . The bulk of the Monte Carlo 
tracks ends up with Lorentz factors 7 > 10 5 , whereas the 
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£ 3 + 
log h [in cm] 



Fig. 1. Electron Lorentz factor as a 
function of altitude for the M74 elec- 
tric field (Eq||) above the surface of 
a neutron star with B = 15.8 x 10 12 
G, P = 0.15 s and R th = 10 5 cm for 
four temperatures of thermal polar 
cap: T 6 = 3.0, 3.5, 4.0, and 4.5. Each 
panel represents 100 curves calcu- 
lated for individual electrons in the 
Monte Carlo way (grey). The mean 
behaviour determined separately by 
solving numerically the differential 
equation (|lj) is overplotted on them 
as the thick solid line. A strong dis- 
agreement between the two methods 
stands out for T§ — 3.5. The dotted 
line is for the case with no radiative 
losses. The thick dashed line sur- 
rounds a region with |7 ICS | > 7 aC c- 
The upper boundary of acceleration 
region for model of Michel (1974) is 
marked by the vertical dashed line. 



integration-method solution gives 7 ~ 64. To understand 
these differences it is worth investigating closely the be- 
haviour of one exemplary electron as determined by the 
two methods. 

In the integration method, the electron's energy in- 
creases monotonically until an equilibrium settles between 
the energy loss process and the acceleration. Fig. |^b shows 
iTics I and 7acc as a function of the Lorentz factor 7 for 
three different altitudes in the case Tq = 3.0. One can 
see that |7 IC s exceeds 7 acc between the two equilibrium 
Lorentz factors 7 oq i and 7 eq 2 the values of which are de- 
termined by the condition |7i CS | = 7 acc - We have calcu- 
lated them using the approximation 7 ICS ~ j Tes where 
7r CS is the energy loss rate due to the resonant part of 
the cross section alone (cf. Eq. (22) in Sturner 1995[ ). The 
two equilibrium Lorentz factors 7 eq i and 7 cq 2 as calcu- 
lated for different altitudes h are presented in Fig. [I] by 
the thick dashed line. For increasing altitude their val- 
ues approach each other because a decrease in the den- 
sity of black-body photons moves down the whole curve 
ITiosCt)! i n Fig. ||b. Thus, the energy-altitude space is di- 
vided into two regimes: the region where |^y ICS | > 7 aC c, 
hereafter called the "energy loss dominated" region, (it is 
surrounded by the thick dashed line in Fig. fjj), and the 



"acceleration dominated" region where 7i C s < 7acc- As 
can be seen in Fig. |l|, according to the integration method 
the trajectory of an electron in the energy-altitude space 
does not penetrate the energy loss dominated region. As 
the electron moves upwards, its Lorentz factor settles at 
the lower equilibrium value 7 oq i, and stays there as long 
as the equilibrium can exist ie. until the curves (7(7) | and 
7acc(7) in Fig. ^b disconnect. For T 6 = 3.0 this happens 
near the altitude h eq = 2 x 10 4 . Above the height h cq the 
acceleration cannot be counterbalanced by the ICS energy 
losses at any Lorentz factor accessible for the electron and 
the increase in the electron energy resumes. 

In the Monte Carlo method, the electron's energy 
losses due to the scattering occur in a discontinuous way, 
thus, 7 oscillates around the value 7 eq i (see Fig. ||). When 
the energy loss rate 7i C s| hardly exceeds the acceleration 
rate (see Fig. |^) the difference between the equilibrium 
Lorentz factors (7 oq i — 7cq2) is small and the electron is 
able to pass through the energy-loss dominated region. 



Eg. for the case T 6 = 3.0, 7 cq 2 



2200 



4x 7 , 



eq 1 



at 



h = 10 3 cm. To increase its energy from 7 oq i to 7 eq 2 it is 
enough for the electron to avoid a scattering over a 'break- 
through'-distance Abt = 7c g, 2 ~ 7cql = 228 cm which is only 

° UL eii|| /mc z 

'a few' times larger than the local mean free path for scat- 
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Fig. 2. Electron energy loss rate due to the ICS (thick 
solid line) and the acceleration rate (thin solid) as a func- 
tion of the electron Lorentz factor for the case Tq = 3.0: a 
case h — with contributions from angular, nonresonant, 
and resonant part of the cross section shown explicitly 
(dotted, dashed, and dot-dashed line, respectively); b de- 
crease in the loss rate with increasing altitude. The curves 
are for h = 0, 10 4 , and 2 x 10 4 cm (from top to bottom). 
Note the narrowness of the Lorentz factor range within 
which |7 ICS | exceeds j acc . 



Fig. 3. Lorentz factor of an electron (a) and the corre- 
sponding local mean free path for the magnetic Comp- 
ton scattering (b) as a function of altitude for the case 
Tq = 3.0. The erratic line is an example of the Monte 
Carlo trajectory. The thick solid line represents a solution 
of the differential equation ()l2|). Note a large range of 7 
for the Monte Carlo track between h = 10 2 and 3 x 10 3 



tering A ICS . The local value of A ICS depends strongly on the 
electron Lorentz factor and for 7 between 7 cq i and 7 eq 2 
it ranges from 28 to 108 cm (the case Tq = 3.0, h = 10 3 
cm). This gives Abt — (8 — 2) x A ICS . Thus, there is a very 
large probability for the electron to gain 7 p> 7 cq 2 at an 
altitude h MC lower than h oq . Once the electron enters the 
acceleration dominated region at h MC <C /i cq , its energy 
starts to increase up to a value much larger than obtained 
in the integration method. 

For increasing temperatures T, the width (7 cq 2 — 7eqi) 
of the energy-loss dominated region increases (Fig. m) 
whereas the mean free path for 7 = 7 cq i decreases. This 
makes the diffusion of electrons through the energy loss 
dominated region more difficult: to increase its energy 
from 7 cq i to 7 Gq 2 an electron must avoid a scattering over a 



distance which is increasing multiplicity of the local mean 
free path. As a result, the energy distribution for outgoing 
electrons becomes softer (Fig. |]). 

It should be emphasized that the strong disagreement 
between the final electron energies as determined with the 
two methods only appears if conditions similar to those 
for the case T 6 = 3.5 in Fig. [I] are fulfilled. These include 
the equilibrium between the maximum rate of resonant 
energy losses and the rate of acceleration at h = /i aC c- 
Making use of Dermer's approximation for |7 res | (Dernier 



1990) one can easily find that it has the maximum at 
7res = e B /(w(l - /3^ min )T) where ^ min = h/ (h 2 + R 2 h ) 1/2 
is a cosine of an angle at which the thermal cap radius is 
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Fig. 4. Energy distributions of electrons at h = 10 cm 
for Tq = 3.0, 3.5, 4.0, and 4.5. The energy unit is mc 2 and 
the convention analogous to F v is used. Both the total 
power of a single distribution and its hardness decrease 
with increasing T. 



seen from the position of electron and w — — In 0.5. The 



equilibrium |7 re 



T 6 cc i = 4.14 



7a 



holds for 



eEn 



1/2 



(1-^Mmin) B 12 {h) 



(13) 



where B 12 (h) = B(h)/ (1Q 12 G) and ^ is in cm" 1 (cf 
Eq. (29) in Chang 1995| ) . For the considered accelerating- 
field (Eq.|l|) this condition becomes 



T 6 eq = 3.46 (3 (1 



-1/2 



(B 12 P) 



-1/4 



(14) 



and is shown in Fig. || for h = and P — 0.15 s as the 
solid line. 

Other conditions required for the diffusion are low 
thermal cap temperatures and high magnetic field 
strengths. Both requirements partially stem from the de- 
pendence A ICS (7ros) oc BT~ 3 which results in A IC s *C Abt 
for high T and low B pc . [Note that at the resonance A ICS 
is proportional to B just as |7 res |]- 

Additionally, high temperatures and low B-field 
strengths preclude the diffusion because of relatively large 
energy losses due to scatterings in the Klein-Nishina 
regime (see Fig 2a). First, for 7 > 7 ros , |7 KN | < 7 acc must 
hold, which requires low T. Second, the resonant bump 
in the I'Vios Ct) I curve must be present, which is possible 
when |7 res (7res)| > |7kn|- Since 7 res (7re S )/7KN oc B, strong 
magnetic fields are preferred. 

For P ~ 0.15 s and Rth > h &cc all these constraints 
limit the diffusion regime to T <; 5 x 10 6 K and B pc ;> 



Fig. 5. Approximate region of strong discrepancies be- 
tween results of the Monte Carlo and the integration 
method for the M74 electric field case (Eq.|l|), for pulsars 
with P ~ 0.15 s and i?th > r pc (hatched). The solid line 
marks the equilibrium between the value of the energy loss 
rate |7ros| near the resonance and the acceleration rate 
7 acc at the surface level (Eq.|lJ with h = 0). The dashed 
line indicates the equilibrium at h — r pc for i? t h = 10 cm 
and the dotted line is for h — r pc , i? t h = 10 4 cm. Dots 
are positions of four pulsars with putative thermal X-ray 
emission (Ogelman 1995[ ). The temperature of a thermal 
cap is not known for B1509— 58: triangles mark the range 
of T considered in Fig. [IJ 



; pc roughly fulfilling Eq. (|14| ) 
taken for h = /i aC c- The region of strong discrepancies 
between the two methods is shown in Fig. || (hatched) . 

Among a few pulsars exhibiting a two-component ther- 
mal X-ray spectra only the Vela pulsar has an inferred 
temperature of the thermal polar cap (T ~ 15.8 x 10 6 K, 



Ogelman 1995) which is sufficient for the acceleration to 
be initially halted by the resonant ICS (ie. it exceeds Tg q 
for h < i?th, see Fig. ||). The equilibrium |7 rcs | = 7 acc 
holds up to the altitude h cq ~ 2.2 x 10 4 cm which is two 
times larger than the inferred thermal cap radius (10 cm, 



Stumer 1995) and two times smaller than r pc = 4.8 x 10 4 



cm. Nevertheless, because of the high T, at h h cq the 
diffusion is precluded by the extremely small A ICS ~ 0.1 
cm. It becomes efficient only at h ;> 0.9h cq , thus, most 
Monte-Carlo electrons follow closely the behaviour pre- 
dicted by the continuous approach by Chang ( |1995| ) or 
Sturner (gggSl). 



For other pulsars with a likely thermal X-ray emis- 
sion from heated polar caps (eg. Geminga, B1055— 52, 
and B0656+14) inferred temperatures are much lower (be- 
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tween 2 and 4 xlO 6 K). Therefore, 7i C s "C 7acc for any 
Lorentz factor accessible for an electron (Fig. ||) and the 
Compton scattering cannot considerably redistribute an 
initial electron energy spectrum. 

Both, in the case of negligible ICS losses, ie. for T <c 
T eq (/i = 0), and in the cases when 7rcs(7rcs)| ^S> 7 acc holds 
up to h — h acc , ie. for T 3> T cq (h = h acc ), we find the en- 
ergy distribution of electrons to be quasi-monoenergetic 
with energy well approximated by Eqjl^. In the case of 
negligible ICS cooling the energy distribution of electrons 
is a narrow peak at the energy j max mc 2 = (E\\h acc , 
whereas for very large efficiency of ICS losses (close to 
100% of j max mc 2 ) nearly all particles are cooled down 
to the energy 7 min for which the ICS loss rate decreases 
sharply (see Fig. ||a and Fig. [I]). Between these limiting 
cases, ie. when the ICS cooling is comparable to the accel- 
eration, broad energy distributions of electrons (ranging 
from 7, nin to 7 max ) emerge (Fig. g). 

3.2. The case of the HM98 electric field 

To enable electron energy losses due to the resonant ICS 
to compete with acceleration (the resonant ICS damping 
becomes then important and the diffusion effect does oc- 
cur) the strength of the electric field E\\ should not ex- 
ceed a critical value which may be roughly estimated as 



3 • 10 4 Bi 2 T|Vcm- 1 (cf. Chang |1995| ). For E\\{h) 



given by eq. (18) of HM98 this may only occur in the case 
of small acceleration length h acc <C r pc (see Eq. @), when 
the acceleration height is limited by the ICS-induced pair 
formation front. Self consistent values of h acc are then of 
the order of a fewxlO 3 cm, i.e. much smaller than for a 
CR-induced case (cf. figs. 5 and 6 in HM98). A value of 
ft-acc is considered self consistent if the acceleration of pri- 
mary electrons by the electric field induced with this h acc 
generates pair formation front at h acc . 

We calculated self consistent h acc in an approximate 
way, following HM98, by finding a minimum of a sum of 
acceleration length and a mean free path for one photon 
absorption: 



(15) 




where to get the first term we used En in the limit 
h < hacc, given by Eq. (§): E\\(h) = Nh with N = 
3rJi? ns B pc KCOSx(c(l - eciO-RD^acc- As an energy of 
pair producing photons in the second term we assumed 
that for the resonant scattering: e ~ ^e B . In Eq. ( |l5| ) 
we neglected a contribution from the mean free path 
for the resonant ICS which is of the order of 150(1 — 
/9 Mmin) -1 Bi2Tq 3 cm near the resonance. Since in the 
regime h acc <C r pc the electric field En depends itself on 
fr-acc we had to repeat the calculation of h acc starting from 
some guess value until convergence. For B pc = 15.8 TG, few percent of electrons reach 7 > 10 2 at h 



P = 0.15 s, pcurv = 10 7 cm and x = 0.1 rad we obtained 
/lace — 2 x 10 3 cm. 

Then we calculated the changes of 7 with h for i? t h = 
10 5 cm and for different thermal cap temperatures (as 
listed in the preceding subsection). We found that the 
most notable divergence between the Monte Carlo results 
and the integration method results occurs for Tq — 3.5. As 
can be seen in Fig. |a, at h ~ 10 3 cm most Monte Carlo 
electrons reach Lorentz factor values between 10 3 and 10 4 , 
contrary to 470 as predicted by the integration method. 
The difference disappears only above the accelerator (i.e. 
for h > /iacc) because such altitude is still low enough (i.e. 
h <C -Rth) for the resonant ICS to damp 7 to its final value 
of about ~ 50 (at which it becomes eventually negligible) . 

An effect similar to that shown in Fig. [IJ appears, how- 
ever, when the resonant ICS losses are no longer effective 
above h = h acc . Fig. ||b shows the case with R t h = 4.0 xlO 3 
cm = 2 x h acc and T§ = 3.0. Broad electron energy distri- 
bution (with average energy much exceeding a value from 
the integration method) emerges. 

Similarly as in the weak electric field case, the position 
of the diffusion region in T-B pc plane (for a fixed value of 
Rth) is here determined with a condition analogous to Eq. 
( |l3| ) (but E\\ of Eq.(||) requires a more accurate determi- 
nation of /iacc) along with a requirement for strong B pc 
and low T. 

For stronger electric fields {E\\ > E^ q ), or for larger 
acceleration lengths the resonant ICS has negligible influ- 
ence on final electron energies. In the first case acceler- 
ation greatly dominates the ICS damping; in the second 
case acceleration is stopped only over a negligible, initial 
part of accelerator height. 

4. Conclusions 

We have calculated the efficiency of energy losses for elec- 
trons accelerated over a hot polar cap of a neutron star. As 
cooling mechanisms we have considered the inverse Comp- 
ton scattering and the curvature radiation, (though the 
latter has been negligible for the considered acceleration 
models). 

The electron energy losses due to the ICS have been 
calculated with the Monte Carlo method and compared 
with the integration method (after Chang 1995| ; Sturner 
|1995 ) based on the prescription for the averaged ICS 
cooling (Dermer 1990| ). We confirm general predictions 
of the integration approach (eg. the limitation of electron 
Lorentz factors below ~ 10 2 for high temperatures), nev- 
ertheless, the "stopping acceleration effect" occurs at a 
slightly higher temperature than that one resulting from 
integrating the differential equation (|l2]). 



For the considered pulsar parameters (B pc = 15.8 x 
10 12 G, Rth = 10 5 cm) and the acceleration potential by 



Michel (1974), the Monte Carlo-based value is approxi- 
mately equal to T MC = 4.5 x 10 6 K (at this T only a 
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Fig. 6. Electron Lorentz factor as a function of altitude 
for the HM98 electric field (Eq.g) with h acc = 2 x 10 3 
cm. Solution found with the integration method (thick 
solid line) is overplotted on ten Monte Carlo tracks. Thick 
dotted line is for the case with no radiative damping. The 
dashed vertical marks h = h acc . a case of i? t h = 10 5 cm 
and T 6 = 3.5; b case of i? t h = 4 x 10 3 cm and T 6 = 3.0; 
the other pulsar parameters are the same as in Fig. [l]. 



in comparison with Ti nt = 3.3 x 10 6 as given by the in- 
tegration approach. Although the temperature difference 
is small, the final electron Lorentz factors as determined 
by the two methods may easily differ by a few orders of 
magnitude (Fig. the case T e = 3.5) if the observed tem- 
perature matches this range. Moreover, it should be kept 
in mind that for most objects with hard X-ray emission 
identified putatively as originating from the hot polar cap, 



temperatures lie within the similar range (2 — 4) x 10 6 K 
(6gelman|l99§). 

We find that for high-S and low-T pulsars the integra- 
tion method gives a poor estimation of electron Lorentz 
factors if the energy loss rate due to the resonant ICS 
hardly exceeds the rate of acceleration. In the case of 
Michel's model and for i? t h > ^acc, this occurs if T ~ 
4(Bi 2 P)" 1/4 which may be the case for PSR B1509-58. 
In the cases when the ICS loss rate dominate firmly over 
the acceleration or is negligible the energy distribution 
of outgoing electrons is quasi-monoenergetic around the 
value which is well approximated with the continuous ap- 
proach. 

We find that preventing electrons from achieving large 
Lorentz factors occurs also for accelerating potential by 
Harding & Muslimov (1998) in its nonsaturated version 



with h acc <C r pc . The discrepancy between the Monte 
Carlo and the integration method also occurs within the 
high-_B and low-T regime and is especially pronounced for 

^?th ~ ^-acc- 
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